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Abstract 



A compound tunneling mechanism from one integrable region to another me- 
diated by a delocalized state in an intermediate chaotic region of phase space 
was recently introduced to explain peculiar features of tunneling in certain 
two-dimensional systems. This mechanism is known as chaos-assisted tun- 
neling. We study its consequences for the distribution of the level splittings 
and obtain a general analytical form for this distribution under the assump- 
tion that chaos assisted tunneling is the only operative mechanism. We have 
checked that the analytical form we obtain agrees with splitting distributions 
calculated numerically for a model system in which chaos-assisted tunneling 
is known to be the dominant mechanism. The distribution depends on two 
parameters: The first gives the scale of the splittings and is related to the 
magnitude of the classically forbidden processes, the second gives a measure 
of the efficiency of possible barriers to classical transport which may exist in 
the chaotic region. If these are weak, this latter parameter is irrelevant; oth- 
erwise it sets an energy scale at which the splitting distribution crosses over 
from one type of behavior to another. The detailed form of the crossover is 
also obtained and found to be in good agreement with numerical results for 
models for chaos-assisted tunneling. 
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I. INTRODUCTION 



Under the denomination of "quantum chaos", a large body of theoretical and experi- 
mental work has been, and continues to be, devoted to the study of the specific features of 
a quantum system which can be traced back to the degree of chaoticity of the underlying 
classical dynamics Among them, questions concerning tunneling effects receive now 

an increasing degree of attention. Indeed, although they are often considered as purely quan- 
tum, since they correspond to classically forbidden events, it appears amply clear now that 
tunneling processes are strongly affected by the nature of the underlying classical dynamics 

This is made most explicit when considering time-dependent systems as well as systems 
in more than one dimension, for which, in contradistinction to the one-d conservative prob- 
lems usually considered in basic quantum mechanics textbooks, energy conservation does 
not constrain the motion to be integrable. A first consequence of having a richer dynam- 
ics is that quasi-degeneracies analogous to those found in the standard symmetric double 
well problem can be observed in systems possessing a discrete symmetry, independently of 
whether a potential barrier actually exists or not. Indeed, unless the dynamics is entirely 
ergodic, some classical trajectories are trapped in d- dimensional invariant manifolds (the 
invariant tori) inside the 2<i-dimensional phase space. If there is a discrete symmetry, say 
parity or time reversal, in a given system, then any tori which are not themselves invariant 
under this symmetry operation will come in symmetrical pairs. Semiclassical EBK quan- 
tization can be applied to these symmetrical tori, and, since they are entirely identical, 
this approximation yields degenerate energy levels. Note that this effect is distinct from 
ordinary symmetry induced eigenvalue degeneracy. Indeed, the effect is only semiclassically 
correct and tunneling effects lift this degeneracy. This kind of tunneling has been dubbed 
"dynamical tunneling " by Davis and Heller [|Kj . 

Rephrasing this in terms of wave functions dynamics makes it clearer why the term 
"tunneling" is proper in spite of the absence of a barrier. Indeed, it is possible to construct 
what Arnold has termed quasi-modes |l l]] whenever there exist invariant tori fulfilling the 
EBK quantization conditions |T^] . These are wave-functions semiclassically constructed on a 
single torus, and which fulfill the Schrddinger equation up to any order in H. Here, similarly 
to what happens when a potential barrier actually exists (e.g. a one-o? conservative system), 
the actual eigenstates are not approximated by a single quasi-mode, but rather by a linear 
combination of quasi-modes constructed on symmetry related tori. Therefore, if one allows 
a quasi-mode constructed on one of the tori, to evolve for a very long time, it will eventually 
evolve into its symmetric partner, whereas classically trajectories remain indefinitely trapped 
on one torus. 

Our interest in higher dimensional systems lies however mainly in the possibility of 
considering dynamics of a different nature. Even for integrable systems, no general theory 
of tunneling in multidimensional systems is presently available. However some theoretical 
and numerical studies clearly demonstrate that the tunneling mechanism in this case 

is rather similar to what is observed for one-d systems. In particular the splitting between 
two quasidegenerate doublets has a smooth exponential dependence in h. 

If on the contrary the dynamics is mixed, (as will generically be the case for low dimen- 
sional Hamiltonian systems), and one is interested in the tunneling between two symmetry 
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related invariant tori separated by a significant chaotic region, new behavior has been ob- 
served, which is quite different from what one is used to in the integrable case. From a 
now growing body of numerical work, either on one- dimensional time dependent systems 
H,|5|||], or on two dimensional conservative systems [7||§, it has become clear that the pres- 
ence of chaos is associated with certain qualitative features, namely: (i) great enhancement 
of the average splitting, (ii) extreme sensitivity to the variation of an external parameter, 
(hi) strong dependence of the tunneling properties on what is going on in the chaotic region 
separating the two tunneling tori. A particularly striking evidence of this last point is given 
in H , where it is observed that reducing drastically the classical transport in the chaotic sea 
from the neighborhood of one torus to the one of its symmetric partners noticeably reduces 
the tunneling rates. 

In U, it has been demonstrated that a natural interpretation of this unusual tunnel- 
ing behavior is obtained if one consider that instead of the "direct" tunneling mechanism 
characteristic of the integrable regime, a new tunneling mechanism takes over, in which the 
chaotic region plays a predominant role. In "direct" tunneling, splittings are caused by the 
overlap of the semiclassical functions constructed on one torus of the symmetrical pair via 
the EBK scheme (the "quasi- modes" ) . As discussed above, these quasi-modes are not eigen- 
states, but fulfill the Schrodinger equation up to arbitrary order in h. Nevertheless, they can 
be connected to other states through "tunneling" matrix elements which are exponentially 
small in h. (See section II. A of || for an explicit example.) In the "chaos-assisted" regime 
on the contrary, the picture is that the particle first tunnels from the integrable region to 
the chaotic sea and from there again to the symmetrical region. More precisely, this means 
that the tunneling is dominated by the matrix elements between the quasi-modes and states 
semiclassically localized in the chaotic region, rather than by the matrix element connecting 
directly the two quasi-modes. This process involves one (or many) intermediate state and 
two tunneling processes. The reason why it may nevertheless dominate is that, since the 
chaotic region is much closer in phase space than the symmetrical integrable region, the 
tunneling amplitudes are expected to be usually much larger. Once in the chaotic region, 
there is nothing to stop the particle from reaching the symmetrical integrable region. In 
this process, not only do the two semiclassical states play a role, but the various delocalized 
chaotic eigenstates which might couple to them also become relevant. For this reason, the 
tunneling amplitudes have a remarkable dependence on h: first, they decay exponentially 
fast with h, reflecting the smooth variation of the tunneling amplitude from the integrable 
region to the chaotic sea. Second, however, and superposed upon this smooth variation, 
there is an extremely irregular fluctuation of the splittings, due to the violent variation in 
the strength of the coupling to the nearby-lying chaotic states. This depends very sensitively 
on the smallness of the energy denominator, i.e. on whether or not a chaotic level lies close 
to the tunneling doublet. In fact, these fluctuations can be so strong as to make any realistic 
assessment of the h dependence of the tunneling amplitude impossible. Note however that 
this picture may be oversimplified: it was observed on the kicked Harper model by Roncaglia 
et al. ||15|| , that direct tunneling can be dominant even in the presence of a significant chaotic 
region. 

That the mechanism described above is actually the one taking place for tunneling in the 
presence of chaos cannot, at the present time, be derived from the basic quantum mechanical 
law of evolution. The numerical evidence as well as a far more careful and detailed discussion 
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of the process are given in || . It is worth stressing, however, that one of the most important 
implications of the chaos-assisted mechanism described above is that it allows for a modeling 
of the splitting distribution in terms of ensembles of random matrices, as discussed in more 
detail in ||. Therefore not only a qualitative interpretation of the numerically observed 
tunneling behavior is obtained, but also quantitative theoretical predictions can be made, 
and compared to numerically obtained data with a surprisingly good accuracy. 

The purpose of this paper being the study of the resulting matrix ensembles, let us be 
more specific on what is understood by "splitting distribution" , and how the theoretical pre- 
dictions have been obtained in ||. As mentioned above, a characteristic feature of tunneling 
in the presence of chaos is the extreme sensitivity of such quantities as splittings, to the vari- 
ation of an external parameter. Within the "chaos-assisted" interpretation of the tunneling 
mechanism, this is quite natural since the splitting of a given doublet may vary by orders of 
magnitudes depending on whether a chaotic state is close to the tunneling doublet or not. 
Therefore, very small changes of external parameters, which leave the classical dynamics 
almost unaltered, may drastically change the splitting by shifting chaotic levels a distance 
of a few mean spacings. Therefore, in any experimental setting the statistical behavior of 
the splittings is likely to be of great relevance, even if one focuses on one single well-defined 
doublet. In this case, the physically relevant quantity will be the distribution of splittings 
for an ensemble obtained by varying an external parameter over a range which is neglible 
on the classical scale, but still large on the quantum scale (i.e. the tunneling doublet crosses 
a lot of chaotic levels). By analogy with the random matrix ensembles describing the spec- 
tral fluctuations of classically chaotic systems^ (see e.g. 0,0), or to their generalization 
introduced to describe partly chaotic systems (see section 5 of 0), one can build random 
matrix ensembles modeling the tunneling distributions corresponding to a given phase space 
structure. 

We shall not repeat here in detail the prescription given in || for constructing the 
random matrix ensemble relevant to a given classical configuration. (A few examples of 
such ensembles are given below.) The following points, however, should be stressed: First, 
we are not considering fully ergodic systems (otherwise there would be no invariant tori); 
hence the chaotic part of phase space cannot a priori be considered as structureless. A 
whole set of partial barriers should quite generally be present, preventing the motion in the 
chaotic region from being completely random. In order to quantify the efficiency of these 
partial barriers, additional time scales must be introduced apart from the mean Lyapunov 
exponent. As a consequence, the random matrix ensembles associated to the chaotic part 
of the phase space cannot be taken as structureless either, as demonstrated in 0. One 
has to introduce "transition ensembles" entirely specified by a set of "transport parameters" 
Ai, A2, • • ■. These transport parameters are fixed by the classical dynamics, and are therefore 
not adjustable model parameters. Therefore the parameters which determine the tunneling 
distribution are, in the chaos-assisted regime, of two kinds: 

— (i) the variance v\ of the tunneling matrix elements. It describes the classically forbidden 



x These were originally introduced in Nuclear Physics without any intention to discuss the nature 
of the underlying classical dynamics 
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part of the tunneling process. There exists, at the present time, no theoretical way to 
evaluate it. 

— (ii) the set of "transport parameters" A 1? A 2 , • • •, which are fixed by the classical dynamics 
inside the chaotic region. 

A second thing we would like to stress is that it is not necessary to solve analytically 
the thus constructed random matrix ensembles to obtain a "theoretical prediction" for the 
splitting distributions. Indeed the latter are entirely specified by the random matrix ensem- 
ble, and can be obtained concretely by performing a rather straightforward Monte Carlo 
calculation: i.e. taking at random a large number of matrices with the distribution specified 
by the ensemble, diagonalizing them numerically to obtain the splitting, and construct in 
this way a histogram of the splitting distribution. This was the procedure used in |§ to 
compare splitting distributions of doublets of a system of coupled quartic oscillators with 
those predicted by the proper random matrix ensemble. 

The obvious drawback of using Monte Carlo simulations to produce the matrix ensem- 
bles prediction is that it sheds little light on the characteristic features of the distribution. 
This is made even more important here because, although all parameters but the tunneling 
amplitude Vt can in principle be computed by studying the classical motion in the chaotic re- 
gion, their practical calculation requires a great deal of effort in the simplest situations, and 
could turn out to be impossible for sufficiently complicated classical structures. Moreover, 
in experimental realizations, the Hamiltonian governing the dynamics may not be known in 
enough detail to allow fixing the parameters of the ensemble with sufficient confidence. 

It is therefore worthwhile to gain some further understanding of the splitting distribution 
determined by the ensembles of random matrices constructed in ||, and to obtain explicit 
expressions for these distributions. As we shall see, they can in fact be expressed in rather 
simple form. Moreover, the resulting distributions have some very specific general feature, 
which can be used as the fingerprint of chaos assisted tunneling even when the precise 
structure of the chaotic motion is unknown. Because of the relative complexity of the 
derivation, we have chosen to organize this paper as follows. In section |I|, we shall give 
the final result, without any justifications, and show how well our analytic findings compare 
with actual splitting distributions, obtained in || for a system of coupled quartic oscillators. 
The remaining of the paper will be devoted to the derivation of this result. Section |T1 
will deal with the simpler case where the chaotic phase space can be taken as structureless. 
In section |V| we shall derive the splitting distributions in the case where effective partial 
barriers are present, and discuss in more details the hypothesis and approximation used in 
the derivation. Section |V| will be devoted to some concluding remarks. 



II. THE SPLITTING DISTRIBUTION 

In the following we introduce some notation. We are interested in a system possessing a 
discrete symmetry P, and for which tunneling takes place between two quasimodes ^ R and 
ty R constructed on symmetrical invariant tori T\ and T2 = P{T\). The eigenstates belong to 
a given symmetry class + or — depending on whether they are symmetric or antisymmetric 
under the action of P. We note and ty R the symmetric and antisymmetric combinations 
of the quasimodes ty R and ty R . If one neglects the direct coupling between the quasimodes, 
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and have the same mean energy Er. The "chaos-assisted" mechanism proposed in || 
assumes that the tunneling from T\ to T 2 originates from the (exponentially small) coupling 



between (resp. ^ R ) and chaotic states of same symmetry \n, +) 



; n 



1,2, 



(resp. 



\n, — )) semiclassically localized in the chaotic region surrounding the islands of stability 
containing T[ and T 2 . Therefore, in a basis where the chaotic part of the Hamiltonian is 
diagonal, the + and - sectors appear respectively as 
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(2.1) 



After diagonalization of H + and H , the regular levels will be shifted respectively from an 
amount 5 + and 5~ , giving the splitting 



6=\5 + - <T 



(2.2) 



Eq. ( |2.1|) merely summarizes the semiclassical picture one has of the chaos assisted tunneling, 
but no random matrix modeling has been introduced yet. This latter is obtained by assuming 
that the statistical properties of the tunneling is correctly reproduced if one model H + and 
H~ by some ensemble of matrices with a specified distribution. 

As discussed in detail in |J, the natural choice for the tunneling matrix elements i>„ 
describing the classically forbidden process is to take them as independent Gaussian variables 
with the same variance vf, and, in the absence of any partial barrier, to use for the chaotic 
levels E^, E 2 , • • • , E±, E 2 , ■ ■ ■ the classical ensembles of Wigner and Dyson which are known 
to model properly the spectral statistics of completely chaotic systems |18j . If time reversal 
invariance symmetry holds, as we shall asssume in the following, this means that one should 
take Ei, E 2 , ■ ■ ■ and E x , E 2 , • ■ ■ as two independent sequences, with a distribution given by 
the Gaussian Orthogonal Ensemble (GOE). Symbolically this ensemble is denoted by 



Er {v} 
{v} {GOEY 



H 



Er {v} 

{v} {GOEY 



(2.3) 



where the subscripts + and — emphasize the independent nature of the distribution. In 
that case, we shall see in section |TTT| that the splitting distribution p(8) is merely a truncated 
Cauchy law 



p(5) 
p(5) 



5 2 + Any? 
(5 > v t ) 



(5 < v t ), 



(2.4) 



As demonstrated in 0, such a simple statistical modeling of the chaotic states does not 
apply any more when structures, such as partial barriers, prevent classical trajectories to flow 
freely from one part of the chaotic phase space to another. In such classical configurations 
(which are presumably generic in systems where chaos and regularity coexist), transition 
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ensembles have to be introduced to obtain a correct statistical description of the chaotic 
states. Compared to the case where no barriers are present, the distribution of chaotic 
states will be modified in two ways. First, each of the parity sequences, taken separately, 
will usually not be distributed as a GOE anymore. However, as stressed in f9j, and as 
will be made extremely clear in section [TTJ, this has a negligible influence on the splitting 
distribution. More important is that such barriers may induce strong correlations between 
the two parity sequences of chaotic states. To fix ideas, let us consider a simple example for 
which a strong partial barrier separates the chaotic sea into two parts IZi and TZ 2 which are 
symmetric images one of each other under P. In that case, the relevant matrix ensemble 
can be symbolically written as || 

H ± -( Er {v} ) (2 5) 

~ \{v} (GOE) s ±(GOE) A (A) ) ' ^ 

where the variance of the matrix elements of (GOE)s is chosen such that it has (in the 
neighborhood of Er) the same mean spacing D as the chaotic states, and the variance a 2 
of the matrix elements of (GOE)a{A.) is fixed by the transport parameter A through 

^ - A . (2.6) 

The transport parameter is in turn semiclassically related to the classical flux $ crossing 
the partial barriers by || 

A = J & (2 7) 

where g — 1/2 is the proportion of states in the corresponding symmetry class, f\ = f 2 = 1/2 
the relative phase space volume of region 1 and 2 and d the number of freedoms. 

For very ineffective barriers, A will be much larger than 1. (GOE)s + (GOE)a(A) and 
(GOE)s — (GOE)a(A) will then be two essentially independent ensembles and one will 
recover in that case the truncated Cauchy law Eq. ( |2.4| ) for the splitting distribution. At 
the opposite extreme, a perfect barrier would give A = (except for classically forbidden 
processes), and the + and — spectra will be strictly identical. In particular so will be the 
displacements 5 + and S~, giving a null splitting. For small, but finite A, a typical level 
E~ will be usually found close to its symmetric analog E+, though slightly shifted. In this 
case, although the displacements S + and 5~ are still distributed as Cauchy, they are strongly 
correlated. This in turn noticeably affects the splitting S. 

Here it may be usuful to give a qualitative description of the meaning of A. It can 
be described as the ratio tu / t c of two timescales, one of which is purely classical, whereas 
the other is essentially quantum mechanical: t c is the time necessary to cross the barrier, 
that is, the typical time that a classical trajectory needs in order to go from one part of 
the chaotic sea to the other. The second time scale tn is the Heisenberg time h/D, where 
D is the average level spacing in the chaotic sea. Eventually, of course, one expects A to 
go to infinity, as the Heisenberg time becomes classically infinite in the semi classical limit. 
Nevertheless there may well be a very large intermediate region, for which A takes very 
small values. In such cases, there exists a classical time scale comparable or larger than 
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the Heisenberg time. This time scale is related to the time necessary to explore all of the 
available phase space. In this respect the situation is quite reminiscent of what happens 
in localization. The difference, of course, is that we only have a small number of weakly 
connected phase space regions and the long classical time has nothing to do with diffusion. 

Usually, the ensemble describing the classical structure of a system will be more compli- 
cated that the simple one given in Eq. ( |2.5|) . Indeed, it is probable that transport from one 
regular island to its symmetric counterpart is affected through the joint effect of a whole 
set of moderately effective barriers rather than by the strong action of a single one. Thus 
one will have to consider much more structured ensembles, with a transport parameter A n 
associated to each barrier. Under these circumstances, there is a possibility that localization 
effects begin to play a non-negligible role. Should localization become as effective in limiting 
the quantum transport as the classically forbidden processes, our treatment would become 



irrelevant. On the other hand we shall see in section [TV] that when this is not the case a 
noticeable simplification of the problem occurs because all the information encoded in the 
transition ensemble (that is, essentially, the A„'s) can be summarized in a single parameter 
a, which is a weighted average of the variance of the (E+ — E~). 

The splitting distribution p(8) therefore depends on three parameters: (ft) 2 , the variance 
of the tunneling matrix elements, a 2 which measures the degree of correlation between the 
odd and even levels, and D, the mean spacing of the chaotic levels to which the are 
connected. If one consider effective barriers, a is smaller than D. Moreover, v t being related 
to classically forbidden processes is usually extremely small, and in particular much smaller 
than a. We shall therefore assume below Vt <C a < D. Then, the main result of this paper 
is that, for this parameter range, the splitting distribution is given by: 

• for Vt < 5 

p(6) = , (2.8) 



for v 2 /a < 5 < v t 

JA1 - 

5 2 + 4irv 2 

for S < vf/a 



P( S ) = 5 2 (Cauchy) , (2.9) 



= 2/rt? , (2.10) 

where the function G is the inverse Fourier transform of exp(— namely 

G{x) = — J exp(iqx) exp(-y/q)dq , (2.11) 

V32r 2 (3/4W 2 

and y.= y- — (2-12) 

7T D z 
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As expected (see section IV. B. of JJ), only the smaller splittings are affected by the transport 
limitation, the distribution for larger splittings being unaffected. The asymptotic behavior 
of the /j,~ 1 G(S//j i ) is given by 



1 D 2 

p(Q) = 2fi- l G(0) = -= 5 in 5 = (2.13) 

I K ! P y J v/2r 2 (3/4) avf K ! 

*~ lG Wri~^^^W- for5 »«^ 2 (2-14) 
Therefore, for small enough a/D and v t /D, the distribution p(S) will, in a log-log plot, 
essentially consist of three straight lines: (i) a horizontal one (at p(0), as given by Eq. ( |2.13p ), 
for < 5 < av 2 /D 2 . (ii) a line of slope (—3/2) in the range av 2 /D 2 < 5 < v 2 /a. (iii) a 
slope (—2) characteristic of the Cauchy distribution for the range v 2 /a < 5 < v t , after which 
the distribution brutally falls to zero. 

Here let us digress shortly to give an intuitive picture of what is happening. Let us first 
assume that there are no barriers. Then chaos-assisted tunneling is essentially a compound 
process involving the two symmetrical states and those chaotic states which lie nearby. Fast 
tunneling (large splittings) occur only if at least one of these states actually lies very near 
to the quasi-degenerate tunneling state. This then yields a tunneling process mediated by 
one single delocalized chaotic state. This process has a characteristic 5~ 2 distribution, as 
we shall show later. Here it is essential to realize that the chaotic state involved, since it 
has a well-defined symmetry, will always directly couple from one torus to its symmetrical 
partner. On the other hand, if we have an efficient barrier, the chaotic states also come 
in quasi degenerate doublets of opposite parity and of width a. Therefore, moderately fast 
processes will be mediated by such a doublet rather than by a single state. Again, we shall 
show that this leads to a universal behavior of 5~ 3 ^ 2 as long as the doublet is identifiable as 
such, that is, as long as the two energy denominators contribute roughly equally. However, 
for very fast processes, the tunneling will again be mediated by a single state, namely the 
one nearest to the tunneling doublet, and the 5~ 2 behavior is again obtained. The details 
are given by the above formulae, which also show a considerable amount of information for 
intermediate cases which cannot be derived in such a simple fashion. 

Before going to the calculation of the above distribution, let us see how it compares to 
actual distributions of splittings obtained numerically for a system of two coupled quartic 
oscillators governed by an Hamiltonian of the form 

2 2 

H(p, q) = P -^~ + a(qt/b + bq 4 2 + 2\ q \ql) . (2.15) 

Except for their presentation (we use here a log-log plot instead of a linear versus log 
binning), the data used in Fig. [I] are exactly the same as those used in Fig. 13 of Ref. 0, to 
which we refer the reader for more precise information on the system investigated. Here, we 
shall only say that each set of data has been obtained by numerically calculating the splittings 
between regular states constructed on a given identified pair of symmetrical invariant tori, 
for various values of the coupling A.@ The range of variation of A is small on the classical scale 



2 In practice, to increase the statistical significance of the distribution, data coming from close 
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(the classical dynamics remains essentially the same), but sufficiently large on the quantum 
scale that a good statistical significance is reached. 

In Fig. |I]we display the comparison between the quartic oscillators data and the predicted 
form of the distribution Eqs. ( |2.8|) - ([2.10|) for two splitting distributions associated to two 
different pairs of symmetric invariant tori. The agreement is extremely good, especially if 
one considers that the distribution extends over more than six decades. Here a remark is in 
order. The parameters a and Vt used for the analytical curves in Fig. [I] are actually tunable 
parameters. This was already the case for v t for the ensemble introduced in |J, since there is 
not yet any semiclassical theory allowing for the calculation of the matrix elements associated 
to such classically forbidden processes. Here however, one has another tunable parameter 
a. In principle, this parameter is fixed once the the random matrix ensemble describing 
the statistical properties of the chaotic level is known, which is the case for this particular 
system. In practice however, there is usually no way to relate a analytically to, say, the set 
of transport parameters A n 's. Therefore a eventually plays the role of a tunable parameter. 
It should be borne in mind however that a and v t only fix the scale of the distribution, 
and in particular the place of the crossover from Cauchy to the G-like behavior, but not its 
shape. Therefore, despite the presence of two tunable parameters, the fact that the splitting 
distribution in Fig. [I] actually follows the prediction of Eqs. (^.8|)-( P7TUD is a very stringent 
test of the relevance of the whole "chaos-assisted" picture. After this attempt to put the 
results in perspective, we turn to their derivation. 



III. THE CASE WITHOUT BARRIERS 



We now want to get down to computing the distribution of the splittings in the case in 
which no barriers are present. Under these circumstances, it is sufficient to compute the 
distribution of 5 + , and hence since these splittings are statistically independent. 

To start with, one should note that, when considering the matrices H + and H~ in 
Eq. ( j2.1|) , the matrix elements vf are associated with classically forbiden processes, and are 
thus extremely small. Therefore, one can compute the displacement 5^ using a first order 
perturbation result. One has however here to take special care with the rare, but important, 
case, where a chaotic level come extremely close to Er. This can be done using the exact two 
by two diagonalization result for each chaotic eigenstate, and adding up the contributions. 
This gives 



N 



R 



Ef 



t 



\ 



\ 



2vf 



Er — Ea 



(3.1) 



(to be understood in the N — > oo limit). In the absence of any values of Ef too close to 
Er, the above equation is equivalent to the usual perturbative result, 



lying tori have been merged. 
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but the full expression Eq. (|3.1| ) has to be used to regularize it whenever any of the (E ± — E R ) 
become of the order of v%. 

Although we shall give below a more detailed discussion of that point, the basic way we 
are going to use Eq. ( |3.1| ) is that the regularized form of 5 ± (Ef,vf) prevents any splitting 
from being significantly larger than vt, and that for 8 ± smaller than Vt, Eq. ( |3.2|) can be used 
safely. To clarify the discussion, we shall for the moment replace the (vf ) 2 by their average 
value t> 2 , and justify below why this does not change the result. Without loss of generality 
we also set E R equal to zero (En is not correlated to the chaotic spectrum, so it can be used 
as the origin of the energies). We shall in addition consider the normalized regular level shift 
x and energy level e\, e2, • ■ • (we drop the superscript + or _ for the normalized quantities) 



x — 2 ; a — p . (3.3) 



With these manipulations, p(x) is in principle obtained for x < D/v t (i.e. 5 ± < v t ) as 
the integral 

p(x) = J 5 (x - P ( e i> e 2> • • • > e Jv) de i de 2 --- de N , (3.4) 

where P(ei, ■ ■ ■ , ejv) is is the joint probability of a GOE spectra with mean density equal 
to one in the center of the semicircle. It appears however that, the correlations of the 
chaotic states have no influence on p(x), because the physics here is determined by the 
singular nature of the energy denominator which is not expected to be very sensitive to 
many-particle correlations. 

To demonstrate this point, let us consider for instance the integral Eq. ( |3.4j) except that 
we take for the chaotic states a Poisson distribution, i.e. that we neglect any correlations 
between them. In that case, introducing £j = iVe*, one can write 

1 N 

iv i=l 

that is the random variable x is the average of the , where the £j's are independent 
variables with density of probability one at the origin. Would the distribution of admit 
a second moment (i.e. (£ -2 ) < oo), the usual central limit theorem would yield a Gaussian 
distribution for p(x). Here, however, the situation is quite different since this variance 
actually diverges. The distribution po(y) of the yi = behaves as y~ 2 for y 3> 1, whatever 
the initial distribution of the £j, as long as that distribution is equal to one for £ equal to zeroS 



3 Therefore, the result will not be affected by a secular change of the mean density of states away 
from the origin. 



11 



From this follows through standard probabilistic arguments [Tj| that x has a non-singular 
limiting distribution, namely the Cauchy law 

= 2 l 2 • ( 3 - 6 ) 
7l 2 + X 1 

Informally, this result can be obtained as follows. If the chaotic states are distributed 
independently, Eq. (|3.4p reads 



roc ^ f 1 N \ 

J -°°i=i \ iV i=i I 



If we now take the Fourier transform of Eq. (|3.7| ), the result factorizes and one obtains 

roo 

p(q) = / p{x)e i<lx dx = p (q/N) N , (3.8) 



where po{q) is the Fourier transform of Po{y)- The large- y behavior of the latter leads to a 
singularity of po(q), which, by reason of the symmetry of Po(y), must be located at the origin. 
Further, this singularity is of the type of a discontinuous derivative. For any symmetric 
function f(y) with the same large y behavior as Po(y), f(y) — Po(y) decreases more rapidly 
that y~ 2 , which implies that f(q) — po(q) has a continuous derivative. The jump in q = 
of the derivative of po(q) must therefore be the same as for f(y) = (1 + x 2 ) -1 , the Fourier 
transform of which is exp(— 7r|g|). Noting moreover that, because of the normalization, 
Po(0) = 1, one has 

Po(q) = l-7r|g| + o(g) (q < 1) , (3.9) 



and therefore 

N 



p(q) = lim (l -J- I = exp(— 7r|g|), (3.10) 



from which the result follows immediately by inverse Fourier transformation. 

At the opposite extreme, one can consider the most rigid spectrum, and see what happens 
if the chaotic states are distributed as a picket fence. In that case, p(x) can be written as 

r+i/2 / +°° -l \ 
p(x)= / deS[x- ( 3 - n ) 



1/2 



which, using the equality 



cotg(vrx) = - + — 53 



X TT X 2 — k 2 

readily gives 

r+l/2 1 

p{x) = dE5{x- 7rcotg(7rE)) = - - , (3.12) 

J -1/2 X z + H z 
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that is the very same Cauchy distribution as for the Poissonian case. There is no doubt 
that if the two extremes of totally uncorrelated and completely correlated spectra give 
the same result, the correlation between chaotic states play little or no role. Therefore, as 
demonstrated in Fig. ||, the splittings are also Cauchy-distributed when the chaotic states are 
GOE distributed. Indeed, this result can actually be shown using supersymmetric techniques 



21| and also turns out to follow from results on ^-matrix ensembles for the one-channel case 
| under quite general conditions for both the Vi and the energies Ei. In fact, it turns out 
that the sums involved in computing the i^-matrix in a one channel system are exactly of the 
type we are interested in and their distribution can be found exactly under the assumption 



that the ensemble is ergodic and analytic in the energy. For details see |22| and |23 |. 

Before turning to the more difficult case of problems where transport limitations play 
a role in the tunneling mechanism, let us come back to a couple of points not treated in 
the above discussion. Since we have seen that correlations between chaotic states are of 
little importance, we shall discuss these points under the assumption that there are no such 
correlations. The first concerns the fact that the tunneling matrix elements are randomly 
distributed following a Gaussian law, instead of being constant as assumed in the above 
discussion. However, it can easily be checked that, in the Poissonian case, this simply 
amounts to performing first the integral over the tunneling matrix elements distribution. 
The second point concerns the need to use the regularized form Eq. ( |3.1| ) instead of its 
nondegenerate approximation Eq. (|3.2|) . Let us now see the effect of using this more correct 
formula which takes quasidegeneracies fully into account by treating the corresponding 2x2 
matrix exactly. In this case the relevant function of & is equal to 



V 




(C N ) 2 - x (CV) 4 +( — I I ; C N = (3.13) 



which is equal to £j~ for & ^> C^ 1 but saturates to a value of l/v t for smaller values of 
This implies that po(q) has a singularity of the type described above only for q less than 
v t /(ND). This in turn involves a departure of p(q) from pure exponential behavior when 
q becomes of the order of Vt/D, and hence for normalized splittings of the order of D/vt, 
which in unnormalized units correspond to splittings of order Vt- We recover in this way 
the intuitive picture discussed above, namely that the basic role of the regularization is to 
prevent having splittings of size larger than the root mean square deviation of the whereas 
the distribution for smaller values remains unaffected. For the distribution p(x) displayed in 
Fig. 0, the effect of this regularization cannot be observed because it only affects the range 
x > = 10 4 which is not covered in this linear scale. It is however clearly seen in Fig. |], 
as well as in the Figs. ^| and [| of the following section. 



IV. THE CASE OF EFFICIENT BARRIERS 

Let us now consider the more complicated case of systems for which transport limitation 
induces strong correlation between the symmetry classes. As we have emphasised in the 
previous section, the correlation between chaotic states inside a symmetry class has little 
or no influence on the distribution of the shifts, 5 + and 5~ , of the regular level, due to 
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their coupling with the chaotic states. Therefore, even in the case where there exist efficient 
barriers to transport in the chaotic region, the shifts S± should be still distributed according 
to the Cauchy distribution derived in the preceding section. In fact, the main effect of 
such barriers is to induce strong correlations between the Ef and the E~, which only affect 
the distribution of the splittings themselves 5 = \5 + — 5~\. Another consequence of the 
presence of partial barriers which may also influence the distribution of splittings is that 
it may yield some inhomogeneity of the variance of the tunneling matrix elements, as well 
as of the correlation between chaotic states. We shall come back to this point at the end 
of the section when comparing our findings with exact results calculated numerically using 
Monte-Carlo techniques. 



A. Derivation of Eqs. (|2T^)-( |2TTo|) 



To lighten somewhat the notation we use in this subsection scaled energies ef = Ef / D 
(for which the mean density of states around the regular level is therefore one), and note 
Vt = Vt/D and at = a/D. We use the following modeling of our problem. First, we shall 
ignore any correlation between chaotic states inside each symmetry class and merely require 
that the mean density of chaotic states be equal to D arround E^ = 0. To normalize the 
number of states to N, we choose the variables ef and e~ both distributed according to a 
Gaussian law, which we take to be N~ 1 e~ 7T( - Ci . In this way, the density at zero is 1/N for 
each level (and thus the total density is 1). Since we shall consider the N —>■ oo limit, this 
can be though of as a flat distribution on the scale of a mean level spacing, the Gaussian 
form being just introduced to normalize in a proper way the distribution. Again, we take the 
chaotic states to be identically distributed random variables, the ef are independent of each 
other for different values of i. The correlations between the two parity sequences of chaotic 
states must, however, be implemented, which we shall do assuming that the (ef — e~) have a 
Gaussian distribution characterized by its width a. As discussed in section [TJ a is related to a 
characteristic time necessary for a classical trajectory to travel from one regular island to its 
symmetric counterpart. To justify this construction, let us consider for instance the ensemble 
of Eq. ( |2.5| ) introduced in section [TT[ In this case, it was seen that the Hamiltonians of the 
two symmetry sectors are related to one another by adding a GOE matrix, the off-diagonal 
elements of which have variance (2a) 2 , related through Eqs. ( |2.6| )-( |277| ) to the classical flux 
$ crossing the partial barrier. As is well-known, the resulting spectrum is formally the 
same as the result of letting the levels moving according to an interacting Brownian motion 
during a time A = a 2 /D 2 , as described by Dyson f24 |. However, for short times (i.e. small 



A), it is generally accepted that an interacting diffusion process can be replaced by a free 
one [p5|l , which here means that the (Ef — E~) follow a Gaussian distribution, of width 
a = \/2(2a). This also follows more specifically from the results shown in p6| . There it 
was shown that a randomly distributed sequence diffusing over a short time A only acquires 
correlations on a scale of vA- (Note that for a 2 ^> D 2 , this modeling becomes essentially 
meaningless. For instance it is not possible anymore to specify unambiguously what E~i is 
to be associated with a given Ef when considering the variance of their difference. This is 
not of great importance, though, since the final result shows that in this case the truncated 



Cauchy behavior described in section III is recovered anyhow.) Thus the joint probability 
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distribution function of the ef and the e„- is 



(^"^(-.(l) 2 ^^-^. (4.1 



P(et,e- 

Note the asymmetric treatment of ef and e~. This simply means that, because of the very 
strong correlation between the two sets of eigenvalues, the large-scale distribution of one set 
entirely determines that of the other. 

For more structured ensembles than the one of Eq. ( j2.5| ), the actual diffusion process, 
while "turning on" the transport parameters A n 's from zero to their actual value, may 
be noticeably more complicated. Nevertheless, because the final splitting results from the 
average effect of the coupling with a large number of chaotic states, it is natural to assume 
that a kind of central limit theorem is involved and that the form Eq. (|4.1| ) can also be 
used in practice (we shall discuss this question in greater detail in the next subsection). As 
mentioned in section |I| however, there will not be anymore a simple relationship between a 
and the transport parameters of these ensembles. 

Now the problem is to compute the splitting distribution. We shall disregard in the 
following the complications created by the inclusion of the complete expression Eq. ( |3.1| ) for 
the shifts, since this only causes a cutoff at values of S equal to v t , as was already discussed 
in the previous Section. As a further simplification we shall take for a moment the tunneling 
matrix elements as being constant (equal to vf) and shall come back later to the (slight) 
modifications to the result due to averaging over their Gaussian distribution. Introducing 
the scaled variable X = (5 + — S~)/D, (S/D = \X\) the distribution of X is obtained as 

N _ 

(4.2) 



p(X) = J ]Q del de~ P(ef, ej) 5 {x 




Again, this integration can be factorized by introducing the Fourier transform F(q) of p{X), 
and everything can be reduced to quadratures. The details are a trifle tedious and are 
therefore relegated to Appendix A. The final result is 

(4.3) 




-q ■ (4.4) 

a 



Here ^(q) is given by the expression 




(4.5) 

poo flf „ 

$( z ) =2 4(1 - cos^e- 1 /* . (4.6) 
Jo t A 



For ease of reference we give the following integrals, which are derived for completeness in 
Appendix B: 
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oo 



*(l/y)dy = vr 3 / 2 (4.7) 
$(l/y 2 )dy = V2^r(3/4) (4.8) 

An asymptotic study of the function for z 3> 1 and 2 < 1 yields 

$(2) = -z 2 \yiz + 0(z 2 ) (*<1) 

= +l + 0(^ 1 ) (z>l), (4.9) 

(the pref actors given here are actually correct, but we shall not need them, and the order 
of magnitude is easy to obtain). Therefore, ^(q) basically gives a measure of the domain of 
y such that q/(l — y 2 ) is larger than one. For q ^> 1, this is obviously of order yfq. It can 
moreover be evaluated more precisely by noting that in this range of q, $ (x/8q/\l — y 2 \^j ~ 

$(V8q/y 2 )- This is the case, because in the range of y where |1 — y 2 \ 1 9^ y 2 , i.e. for y not 
large, both q/\l — y 2 \ and q/y 2 are much larger than one if q is, and therefore $ saturates 
to its asymptotic value 1 anyhow. One therefore finds for q ^> 1 



~ f°° f^F-J = 2 3/4 r(3/4)y / 2^. (4.10) 

j-00 \ y J 



and therefore 



F(q) ~ F^g) = exp (-2 3 / 4 r(3/4) yjav?\q\) for g > -| . (4.11) 

For 9 <C 1 on the other hand, g/(l — y 2 ) is large only in the neigborhood of y = 1. From 
this follows that one can restrict oneself to the range of integration y ~ ±1. One obtains 

fl+q 

ty(q) « 2 / 

Jl—Q 



DC 



2 /" =(2vr) 3 / 2 g, 



00 



and therefore 

F(g) « F (g) = exp (-2™ t 2 |g|) for g « ^; . (4.12) 

It remains to perform the inverse Fourier transform Eq. (|4.4| ), and to deduce the asymptotic 
behavior oip(X) from the one of F(q) . For large X, p(X) is dominated by the singularities of 
its Fourier transform, which here means the derivative discontinuity at the origin. Therefore 
for v t > X ^> v 2 /a one can use in Eq. ( |4.4| ) the asymptotic q <C 0/v 2 approximation 
F (q) of F(q). Applying the inverse Fourier transformation yields an almost perfect Cauchy 
distribution of the form 



= ' (413) 
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The above result amounts to adding independently the variable 5 + and S~, distributed as 
given by Eq. (|3.6| ), and to neglect the correlations between the two symmetry classes. As 
mentioned in |J, this is indeed quite natural since splittings larger than vf/d are due to 
chaotic levels coming closer than a distance at of the regular level. Since at can be viewed 
as the scale on which chaotic levels are correlated, chaotic states contributing to p(X) for 
X ^> vf/a can therefore be considered as essentially decorrelated from their symmetric 
counterparts. In the language of section 2, the splittings we are looking at here are so large 
that the tunneling is always mediated by a single state. 

Let consider now the range X <C v 2 /a, for which the splitting distribution is affected 
by correlation between symmetry classes. In that case, the term exp(— iqX) in the inverse 
Fourier transformation Eq. (|4.3| ) is essentially constant in all the range q < a/v 2 where F(q) 
differs from its asymptotic behavior Foo(g). Noting that 

±- ( r F^q) exp(- iq x)dq) = A^G^/A) ; A = V8T 2 (3/A)av 2 , (4.14) 

2,71 \J-oo / 



where G(x) is the inverse Fourier transform of exp(— y/q) as defined in Eq. ( [2.11| ). One 
therefore has 

p (X) = \- l G{X/\) + K (4.15) 

where K is the constant 

K=±- r dq{F{q)-F 00 (q)) (4.16) 

ZlT J-oo 

For small a however, K is of order a 2 /v 2 , when A _1 G(X/A) range from order l/(av 2 ) at 
X = to a 2 /v 2 at its lowest value, i.e. at the crossover X ~ v 2 /a between the Cauchy- 
like and G-like behavior. Therefore it can usually be neglected, although in some special 
circumstances it shows up as a small plateau between these two regimes; we shall disregard it 
from now on. Then, the large- AT behavior of A _1 G(X/A) is dominated by the y/q singularity 
at the origin of F^q), so that it goes as X~ 3 / 2 as X — > oo. This is in fact the hallmark of 
the X <^ v 2 /a regime we are discussing. 

Finally, one has to take into account the fact that the tunneling matrix elements are not 
constant, but randomly distributed. As can be seen in the derivation of Eq. (|4.4j ) (see the 
remark below Eq. ( |A4| ) ) this merely amounts to replacing the expression for ^(q) given in 
eq. fl4.5|) by the function \I/(g) defined as follows: 

= (A/a) ) , (4-17) 



where the brackets denote averaging over the v 's, which we take to have a Gaussian distri- 
bution with variance v 2 . This new function is of course much more complicated than the 
original one, but its asymptotic behavior for small or large values of q is readily obtained. 
Indeed, for g<C 1, ^(q) is proportional to |g|, so that the average is obtained by replacing v 
by v t . On the other hand, for q 1, ^(q) is proportional to \f\q\, so that its average over 

a Gaussian distribution is obtained by replacing v by ^2/7iv t . Using these facts together 
with the above estimates for the behavior of p(X) one finally obtains the result stated in 
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Eqs. ( |2.8j - |2.1(3|) in section |T[ Some trivial differences: There we consider splittings as being 
always positive, whereas in the above computation we treated positive and negative split- 
tings separately: this introduces a new factor of two. Further, we have made the dependence 
on D explicit, which in particular means replacing p(X) by Dp(X/D) as well as replacing 
a and v t by their original expressions. 



B. More Structured Ensembles 

For the simple ensemble Eq. (|2.5|) , the two main assumptions we made to replace the 
exact distribution by Eq. ( |4.2| ), namely to neglect correlation of chaotic states among a given 
symmetry class and to replace the interacting brownian diffusion by a free one for small a are 
really under control. Indeed, section |T| gives full justification of the first assumption, and the 
second can be seen as a simple consequence of standard perturbation theory. And actually, 
as shown in Fig. |3], one can see that our analytical findings perfectly agree wit an "exact" 
Monte-Carlo evaluation of the splitting distribution generated by the ensemble Eq. ( |2.5| ). 
We stress that in this very simple case the parameter a that we are using is (for small a) 
simply related to the variance a 2 of the non-diagonal matrix elements of (GOE)a (indeed 
a 2 = 2(2cx) 2 ), and that therefore there are no adjustable parameters in this comparison. 

More structured ensembles deserve however some further discussion. Let consider for in- 
stance the ensemble relevant to the quartic oscillator system used as illustration in section ||. 
Symbolically, this ensemble can be written as M 



( E R {v} \ 

{v} (GOE), (GOE)±(A 12 ) (GO J E) ± (A 13 ) 

(GO J E) ± (A 12 ) {GOE)f 

V (GOE)±(A 13 ) (GOE)f J 



(4.18) 



where the subscript ± again indicates ensemble which are independent in the + and — 
symmetry class. Noting D tot the total density of states (in a given symmetry class), (GOE)i 
stands for a Gaussian Orthogonal Ensemble such that the mean level density in the center 
of the semicircle is fiD tot , and (GO£') ± (Aj fc ) represent Gaussian distributed independent 
matrix elements of variance a 2 k = AjkD 2 ot . (For the configuration of the quartic oscillators 
corresponding to Fig. || one has f\ = 0.5, f 2 = 0.2, / 3 = 0.3, and A i2 = 0.14, A i3 = 0.11.) 
For such complicated ensembles the two assumptions concerning the irrelevance of intra-class 
correlations and essentially Gaussian distribution of the (Ef — E~) are presumably equally 
well fulfilled here as in the simple case of Eq. ( |2.5| ) . What is lost however is the uniformity 
of the distribution of the tunneling matrix elements and of the variance of the (Ef — E~). 
Indeed, in the above example, a diagonalization of the chaotic part of the Hamiltonian is 
going to transfer some tunneling matrix elements from the block connecting Er to (GOE)i 
to the ones connecting (GOE) 2 and (GOE) 3 . One may end in this way with three different 
scales for the variance of the tunneling matrix elements as well as for the parameter a (one 
for each (GOE) block). 

More generally, the typical situation will be that a [possibly large] number of (GOE) 
blocks, (GOE)i, (GOE) 2 , (GOE)k are involved in the tunneling process. After diago- 
nalization of the chaotic part of the Hamiltonian, both the variance of the tunneling matrix 
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elements, and the degree of correlation between symmetry classes, will be block dependent. 
Each block (GOE) k (k = 1,...,K) would have then to be characterized by a tunneling 
parameter v k and a transport parameter a k (a k and v k highly correlated), in addition to its 
dimension N k = f k N. Let us introduce the notation 

I(a,v t ;q) = -——V^q/aD) , (4.19) 
JJy Z7T 

where ^(q) is defined by Eqs. ( |4.4| ). A straightforward modification! of the derivation of 
Eq. ( f4.4| ) gives that taking into account the block dependence of a k and v k merely amounts 
to replacing this equation (i.e. F(q) = exp(— I(a,v t ; q))) by 

F{q) = exp (- £ f k I{a k , v k ; q)) . (4.20) 

Inspection of Eqs. (|4.11|) and (|4.12|) then shows that they remain valid provided a and Vt 
are defined now as 

^ = Y.M ( 4 - 21 ) 

k 

a = \(j2f k al /2 v k ) . (4.22) 
v t \ k / 

Multiplying Eq. fl£21p by N, Nv 2 t appears as the [average] square norm of the projection of 
the quasimode on the chaotic space. It is therefore independent of the chaotic phase space 
structure. This, for instance, allows to compute v t from the variance of the tunneling matrix 
elements before diagonalization of the chaotic part of the Hamiltonian. The parameter a 
and Vt have moreover a certain number of intuitively clear properties: If all v k are multiplied 
by a constant factor, the effective tunneling element Vt is multiplied by the same factor, 
whereas a is unaffected. Further, if all v k are identical, then the effective tunneling element 
is the same. On the other hand, the same is not true of a: If all a k and all v k are taken 
to be equal, the effective efficiency of the classical barrier now depends on the number of 
different components of phase space through which tunneling can take place. Further, we 
see that any components with negligible values of v k will contribute negligibly both to v t 
and a. Thus we can identify a given part of phase space through which tunneling actually 
occurs and limit ourselves to it. 

With the definitions Eqs. ( |4.21|) and (|4.22 ) of Vt and a, structured ensemble are therefore 



seen to behave in essentially the same way as the simple ensemble Eq. (|2.5|) . The, only 
difference is that the condition of validity of Eqs. (|4.11|) and fl4.12|), that is respectively 



In Eq. (CO]), lmiAr-^ooCl - Hq)/N) n has to be replaced by 



lim 



K 



H(l-I(a k ,v k ;q)/N) 
Lk=l 



fkN 
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and 



q > (a k D)/vl 



(4.23) 
(4.24) 



must now be fulfilled for all fc's. The transition between the different regimes of the distri- 
butions may therefore be less sharp than for the ensemble Eq. ( |2.5|) . 

If the partial barriers structures were to become highly developed, say to the point that 
the ensemble could meaningfully be described in terms of band matrices, then obviously the 
issue of localization should have to be considered. In this case, the orders of magnitude of 
the «fc might become comparable to those of the Vk, and most of the splitting distribution 
might be actually in a transition-like regime. This is exactly the sort of problems we pointed 
out in our earlier discussion of the physical situation. However, as long as the a^'s are clearly 
larger than the i^'s, the transition from one regime to another should still take place on 
a short scale as compared to the range spanned by the distribution. Physically speaking, 
this condition amounts to saying that classically forbidden processes are always much slower 
than classically allowed ones. In that case, the form of the result should not be noticeably 
affected. For example, as seen is Fig|], the distribution resulting from the ensemble described 
in Eq. (4.18) still perfectly follows the predicted form Eqs. ( |2.8| )- Q2"TT0| ) (note however that 
a is now a tunable parameter). 



As a conclusion, we have provided in this paper an analytical study of the splitting 
distributions generated by ensembles of random matrices constructed in [[| to model a 
tunneling process in the chaos assisted regime. The original ensembles may contain such a 
complicated structure that a general answer to this problem may seem a priori out of reach. 
Nevertheless, it turns out that only the average size of the tunneling matrix elements and 
the degree of correlation between the chaotic spectra in the different symmetry class affect 
the distribution, and that therefore the problem can be reduced to a simpler formulation 
which is tractable. 

The basic reason for the considerable simplifications encountered was in essence already 
pointed out in ||. It is due to the fact that for large splittings only the situation of near- 
resonance to a given state of the chaotic sea is of relevance. To this obvious remark, we only 
need add that for the case in which efficient barriers are at work, the tunneling operates not 
through single states, but through quasi-degenerate doublets of states of opposite parity. 
These are of course less efficient in promoting tunneling, since the particle requires a time 
of the order of the width of the doublet to reach the symmetrical torus. In either case, the 
behavior is determined by rather natural probabilistic considerations. It turns out to be 
sufficient to consider only the probability of one single eigenvalue being near the tunneling 
state, so that correlations between eigenvalues and the like could be safely ignored. Further, 
the very simplicity of the physical picture given here results in it being fairly robust to 
changes in minor details of the model. Thus it does not appear necessary that all states in 
the chaotic sea should participate equally in the tunneling process, nor that the couplings 
should be uniform. In fact, the main limitations of our result seem to be the ones related 
to localization phenomena. If the structure of the barriers in phase space is sufficiently 
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complicated, it is possible that localization effects, associated to the presence of a large 
number of partial barriers, become as effective in limiting tunneling from one quasimode to 
its symmetric partner as the initial classically forbidden process. In that case, the splitting 
distributions we have obtained would not be relevant anymore. However, this should not 
be a too severe limitation, and it should generally be possible to determine for any given 
system whether this takes place or not. When it does not, the picture of tunneling in the 
presence as well as in the absence of barriers to transport is indeed the one we gave. This 
is substantiated by the numerical work done: In particular, we showed that not only the 
simplest model of a barrier gives results in good agreement with theoretical predictions, but 
also a highly specific random matrix ensemble constructed explicitly in order to model chaos- 
assited tunneling in a system of coupled quartic oscillators was well fitted by the theoretical 
predictions, as were also the actual splitting distribution for this system. 

This might possibly open up a way to identify chaos-assisted tunneling in experimental 
systems. In such systems, the exhaustive study of the classical mechanics necessary to 
produce a satisfactory random matrix ensemble would probably not be feasible. Nevertheless 
the above remarks strongly suggest that if chaos- assisted tunneling is present, the splitting 
distribution will reflect the fact by showing a highly specific and well-characterized behavior. 
Indeed, as discussed throughout this paper, only the scale of the distribution and the position 
of the transitions between the different regimes are system dependent, but the shape of the 
distribution is essentially universal. In particular the experimental detection of a transition 
from a S~ 3 ^ 2 behavior, characteristic of the G-like regime, to a S~ 2 behavior characteristic 
of the Cauchy-like regime, would be a powerful argument in favor of the presence of chaos- 
assisted tunneling. 



We acknowledge many helpful discussions with Oriol Bohigas, as well as with Steven 
Tomsovic which we thank in addition for providing us with the numerical data of the quartic 
oscillators system used in Fig. |l|. Denis Ullmo wants to thank the Institute for Nuclear 
Theory, Seattle WA, where parts of this work has been done. Frangois Leyvraz would like to 
thank the Division de Physique Theorique for its kind hospitality during his sabbtical stay 
there where most of this work was done, as well as acknowledge financial support by DGAPA 
during that same period. The Division de Physique Theorique is "Unite de Recherche des 
Universites Paris 11 et Paris 6 Associee au CNRS". 

APPENDIX A: COMPUTATION OF THE DISTRIBUTION FUNCTION 



Denote by brackets the integration over ef and e i with the weight function P{e'l,e i ) 
(see Eq Q. We define 
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F(q) = / F(X)e iqX dx 
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'expLv^i^-^))) (A2) 



This last expression factorizes in A?" factors, each of which is a double integral. Denoting the 
corresponding average over e + and e~ also by brackets, one obtains: 



F(q) = (exp (iqv 2 t {l/e + - 1/e"))) (A3) 
= (l - (l - exp (V) t 2 (l/e + - l/O)))^ 



N 



N 



where the last line defines In{q)- The reason for this manipulation is that in this way In 
goes to a finite limit I as N —>■ oo and therefore 

lim F(q) = lim (1 - I(q)/N) N = exp(-I(q)) . (A4) 

Note moreover that taking into account the fact that the tunneling matrix elements are 
random variable of variance v 2 instead of being constant just amount to understand (■) as 
containing a further integral on the tunneling matrix elements distribution. This introduce 
no further difficulties in the calculation of J, except for still heavier notations. We shall 
therefore not consider it in this appendix, and just modify the final result in the appropriate 
way at the end of section [IV A . 
One finds 

I{q) = -jL- J de + de- (l - exp («/*g(l/e + - 1/e-))) exp (- {e+ ~J H , (A5) 



since, the above integral being convergent, liniTv^oo In is just obtained by dropping the term 
—7i(e + /N) 2 in the exponent of P(e + , e~) (see Eq. ( f4.1| )). Making the successive transforma- 
tions y = (l/e+ + l/e-)/(l/e+-l/e-), w = (l/e + -l/e-), followed by t = w(a(l -u 2 )/V8), 
I(q) can be expressed as 



/W = ^/7^/^««-l)e X p 



2ira J (1 — y 2 ) 2 J |w| 3 \ a 2 (l—y 2 ) 2 w 2 



1 r dyfg-(l-cos(^Q)e-^. (A6) 



y/2n J J \t\ 3 V \a(l-y 2 ) 
If we now introduce $(z) as in the text, 

poo pjf , „ 

$M = 2 ^(1 - coazQe-V* , (A7) 

Jo t 6 

one easily obtains 

at f j ^ ( \^8v 2 q x 



^ = wJ iy *W^\l (A8) 
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From this follows the formula given in the text: 



F(q) = exp f - -jL *(v*q/a) J , (A9) 



where ^(q) is given by the expression 

* te K>*(h^k) • < A10 > 



APPENDIX B: SOME USEFUL INTEGRALS 

We first give another expression for $(?/): 

$(y) = l°°dt(l- cos-^e"*, (Bl) 

which is obtained from the original definition by substituting 1/t 2 by t. From this follows 

dy 



r$(i/ y )dy= r m% 

J-oo J-oo y z 



poo POO 

= / dte~ f / (1 - cos 

JO J-oo 



00 dt . f°° 1 - cos y , 

e"* / (B2) 



^-00 y 2 

= r(i/ 2) r *y!*=^> 

J-00 V 



-00 y 

The other integral is handled similarly: 

dy 



n<S>(l/y 2 )dy= f°° ®{y 2 )% 

J —00 J —00 v 



00 dt _ t r°° 1 - cos y 2 
J-00 f 



e-< ^^dy (B3) 



/oo 
siny 2 A/ = v / 27r(3/4). 
-00 
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FIGURES 



FIG. 1. Comparison between the quartic oscillator's tunneling splitting distribution (square 
symbols) and the predicted form Eqs.( |2.8D -( |2.10 ) for two different groups of tunneling tori. Except 
for their presentation, the quartic oscillator's data are the same as those in Fig. 13 of Ref. @. The 
transition from the G-like behavior (solid line) to Cauchy-like behavior (chained dot) is clearly 
seen, in spite of this latter being valid on a much shorter range, (a) Group To (using the notations 
of Ref. with vt = 1.1 10~ 2 . (b) Group T\ , with vt = 2.5 10 -2 . The transport parameter has 
the same value a = .04 in both cases, consistent with the fact that the partial barriers structure 
is the same in both cases. It has been taken into account that only a fraction D c r = 0.36Z) of the 
states are actually participating to the tunneling process. 



FIG. 2. Comparison between a Monte Carlo calculated distribution of the reduced variable 
x = 5 ± (D/v 2 ) (solid line) and the Cauchy law Eq. Q3.6| ) (chained dot). The Monte Carlo result is 
obtained from the numerical diagonalization of 10 5 matrices of size 80 x 80, which matrix elements 
are taken at random with the distribution specified by the ensemble Eq. ( |2.3| ) (using vt/D = 10~ 4 ). 
It thus take fully into account the GOE correlations of the chaotic states. Nevertheless, and 
although a linear scale has been used to emphasize the center of the distribution where the effects 
of correlations should be the strongest, the two curves are essentially undistinguishable. 



FIG. 3. Comparison between a Monte Carlo calculated distribution of splittings 5 for the 
simple ensemble Eq. ( |2.5| ) (solid line) and the the predicted form Eqs. fl2.8|) -( [2.10 ). The parameters 
of the Monte Carlo calculations are A = 10~ 2 /8 (imposing a/D = 0.1 for the theoretical curve), 
vt/D = 10~ 4 , number of matrices: 3.10 5 , size of matrices 60 x 60. The three regimes: G-like 
behavior (chained dot), Cauchy-like behavior (dash), and truncation of the Cauchy law for splitting 
greater than vt/D are clearly seen. 



FIG. 4. Comparison between a Monte Carlo calculated distribution of splittings 5 for the 
ensemble Eq. ( PI) with /i = 0.5, / 2 = 0.2, / 3 = 0.3, and A i2 = 0.14, A i3 = 0.11 (solid line), and 
the the predicted form Eqs.( |2.8] )-( 2.1C| ). The Monte Carlo calculations has been performed with 
10 5 matrices of size 100 x 100, using as tunneling parameter (vi) 2 /D = 10 -3 . For the theoretical 
curves, namely the G-like (chained dot) and Cauchy-like (dash) behaviors, the tunneling parameter 
is determined by Eq. ( 4.21[ ) as (vt) 2 = (v\) 2 /2. The transport parameters is however here a tunable 
parameter, which has been taken equal to a/D = 0.1. 
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